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ABSTRACT 



Context. 

Aims. We search for the most suitable set of cosmological parameters that describes the observable universe. The search includes the 
possibility of quintessential flat universes, i.e., the analysis is restricted to the determination of the dimensionless matter density and 
the quintessential parameters, Qu and wq, respectively. 

Methods. Our study is focused on comparing the position of features at large scales in the density fluctuation field at diff'erent redshifts 
by analysing the evolution of the quasar two-point correlation function. We trace the density field fluctuations at large scales using a 
large and homogeneous sample of quasars (~ 38000 objects with 0.3 < z < 2.4 and a median z - 1.45) drawn from the Sloan Digital 
Sky Survey Data Release Six. The analysis relies on the assumption that, in the linear regime, the length scale of a particular feature 
should remain fixed at diff'erent times of the universe for the proper cosmological model. Our study does not assume any particular 
comoving length scale at which a feature should be found, but intends to perform a comparison for a wide range of scales instead. 
This is done by quantifying the amount of overlap among the quasar correlation functions at diff'erent times using a cross-correlation 
technique. 

Results. The most likely cosmological model is CIm = 0.21 + 0.02 and wq - -0.93 + 0.04, in agreement with previous studies. These 
constraints are the result of a good overall agreement of the correlation function at diff'erent redshifts over scales ~ 100 - 300/z~^Mpc. 
Conclusions. Under the assumption of a flat cosmological model, our results indicate that we are living in a low density universe with 
a quintessential parameter greater than the one corresponding to a cosmological constant. This work also demonstrates that a large 
homogeneous quasar sample can be used to tighten the constraints upon cosmological parameters. 

Key words, cosmological parameters - cosmology: observations - large scale structure of the universe - quasars: general 



1. Introduction 

Observations of the cosmic microwave background (CMB, e.g. 
Spergel et al.l |2Q07 | ) and type la supernovae (e.g. iRiess et all 
1998l:lPerlmutter et al.|[l999l:lAstier et al.ll2006h support the idea 
that the missing energy in the universe should possess negative 
pressure p and an equation of state p = wp. One possible can- 
didate for the missing energy is the vacuum e nergy density o r 
cosmological constant A for which w = -1 (iPeeblesI [l984l) . 
The resulting cosmological model, ACDM, consists of a mix- 
ture of vacuum energy and cold dark matter. Another possibility 
is the QCDM cosmolo gy based on a mixture of quintessence 
and cold dark matter (Ra tra & Peebles|[l988t) . The quintessence, 
which is defined as a fifth element diflTerent from baryons, neu- 
trinos, dark matter and radi ation, is a slowly -varying spatially in- 
homogeneous component (ICaldwell et al .1 11998) whose w value 
is less than 0. Many studies restrict the range of w to the inter- 
val (-1 < w < 0) since this range best fits current cosmological 
observations. 

The distinction between ACDM and QCDM models is 
of key importance in cosmology. Even when both models 
match the CMB observations well, the QCDM models seem 
to be more suitable to describe the universe at high redshifts 
(lEisenstein e t al. 2005). According to deep redshift surveys, the- 
ories should predict a strong large scale clustering structure 
and quasar formation at this stage. For this to happen, the 



cessation of growth of the structures should occur at earlier 
times. Unlike ACDM, QCDM universes satisfy this requirement 
since the larger the values of w, the earlier the growth ceases. 
Nevertheless, it is necessary to know how big the diflTerences 
among both cosmologies are in order to restrict the possible ini- 
tial conditions in the observed universe. To do so, a tuning pro- 
cess of the cosmological parameters is needed, combining the 
best current observations with the most suitable statistical tools. 

Several attempts have been made to tighten the constraints 
on cosmological parameters. Among the tests suggested in 
the literature w e can mention measurements of the Hubbl e 
parameter (e.g., IWei & Zhang j2QQ8: Szvdtowsk ietaD l2QQ8h 
gravit ational lens ing (p - g-. IZhang et al.. .2007: iZhu & Serend 
120081: iDore et al.l l2007l: iLa Vacca & Colombo! [2008h. angular 
sizes of distant objects (e.g.. iDaly et al.l 120071: Santos & Lima 
20081) , gas mass f raction in galaxy clusters (e.g., [Ch en & Ratrd 



20041: ISenI l2008l) and baryon acoustic oscillation peaks at 



lar ge scales (e.g., iLima et al .1120071: ISapone & Amendolal[2QQ7l: 
San chez et a"l]l20Q8h . 

In this work, we rely on the fact that if a particular spa- 
tial scale can be measured at diflTerent stages of the evolution 
of the universe, that scale can then be used to determine the 
values of the cosmological parameters. At large enough scales 
(r > 100/z"^Mpc), density perturbations are in the linear regime 
and above the present-day turnaround scale. Thus, structures at 
large scales, for instance features in the two-point correlation 



2 



Zandivarez & Martmez: Constraining flat quintessential universes 



function, should be fixed in comoving length scales at diff'er- 
ent times in the evolution of the universe. Since quasars can 
be observed at very large distances, studying the evolution of 
quasar clustering, i.e., identifying similar features in the density 
field at diff'erent redshifts, could restrict the possible universes 
that match the observa tions. This approach has been adopted by 
iRoukema et al.l (l2002h , hereafter ^02, using a relatively small 
sample of quasars and finding a wide range of possibilities for 
the quintessential parameter (-1.0 < w < -0.5). 

We use a large sample of quasars drawn from the Sixth Data 
Release (DR6; Adelman-McCarthv et al. 2008) of the Sloan 
Digital Sky Survey fSDSSi lYork et al.ll2000i) to infer cosmolog- 
ical parameters. Our study relies on the analysis of features at 
large scales in the quasar correlation function. Since there is a 
large amount of evidence in the literature that we are living in a 
flat universe (e.g. Sanchez et al. 2006), we are interested in con- 
straining perfectly flat quintessential universes, i.e., in our anal- 
ysis there are only two parameters to be determined: the dimen- 
sionless matter density parameter flu, and wq (w denoted with 
the subscript Q to specify QCDM cosmologies). It is important 
to note that QCDM cosmologies allow a wide spectrum of pos- 
sibilities for the equation of state wq, which could be constant, 
uniformly evolving or oscillatory. In this work, we restrict our 
analysis to QCDM universes where wq is constant since there 
is evidence that this parameter does not evolve in time (e.g., 
IWang & Muk heriee 2007; Gong et al. 20M)- We follow the line 
of thought of R02 in which they did not assume the scale at 
which these features should occur, but only required consistency 
of the scale between diflTerent times. 

This paper is organised as follows: in Sect.[2]we describe the 
sample of DR6 quasars; in Sect.[3]we constrain the values of the 
cosmological parameters Qm and wq, analysing features in the 
quasar correlation function at diflTerent redshifts; we discuss our 
results in Sect. IH Throughout this paper we assume a Hubble 
constant Hq = 100 h km s"^ Mpc"^ and all magnitudes are in 
the AB system. As stated previously, we also assume that the 
universe is perfectly flat, i.e, Dm + = 1- 



2. The quasar sample 

The SDSS DR6 provides the largest homogeneous sample of 
quasars currently available. The identification of quasars is 
achieved with high efl&ciency and com pleteness thanks t o the 
wide-field five band photometric system (iGunn et al ]fT998h . The 
high homogeneity of the SDSS spectroscopic sample is ac- 
complished by selecting the spectroscopi c targets consistently 
on the basis of th eir photometric data (iRichards et al.l I2OO2I 
iBlanton et alll2003h . 

To ensure further homogeneity in our quasa r sample, we se- 
lect quasars that match the following criteria (ISchneider et al.l 
|2005; Yahata et al. 2005): 

- The quasar primTarget flag is either "QSO_CAP" 
or "QSO_SKIRT" or "QSO_FIRST_CAP" or 
"QSO_FIRST_SKIRT"; 

- /-band PSF mag nitude corrected for G alactic reddening (us- 
ing the maps bv ISchlegel et al.lll998l) satisfies 15.0 < / < 
19.1 which are the limits of the quasar target selection algo- 
rithm; 

- The redshift is z < 2.4, since the incompleteness becomes 
important for higher redshifts; 

- The absolute magnitude is brighter than Mi = -22.775 - 
51og(/z), computed assuming Dm = 0.3, wq = -1.0 and a 
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Fig. 1. Quasar redshift distribution shown as thin line histograms 
in both panels. The upper panel shows in dotted lines the five 
Gaussian functions whose sum (thick line) fits the redshift dis- 
tribution of DR6 quasars. The lower panel shows in thick line 
the redshift distribution of our random quasar sample obtained 
by using the fit shown in the upper panel. The vertical dashed 
line shows the 50th percentile (median value) corresponding to 
our two redshift bin analysis, while vertical dotted lines show 
the 25th (first quartile) and 75th (third quartile) percentiles used 
in the four redshift bins case. 

^-correction corresponding to a power-law quasar spectrum 
with a - -0.5; 

- We also exclude the southern Galactic region of SDSS due 
to diff'erences in the selection compared to the main survey. 

Our final quasar sample comprises 38060 quasars with a median 
redshift Zmed = 1-45. The redshift distribution of these objects 
can be seen as a solid line histogram in Fig.[Tl 

3. Determining tlie parameters Qm and wq 

The key point of our work is to search for features in the redshift- 
space two-point correlation function of quasars that remain fixed 
in comoving scales as the universe evolves. We do not attempt to 
model theoretical quasar correlation functions to compare with 
observations since this would require taking into account sev- 
eral not- well-known factors such as redshift space distortions, 
quasar bias, survey selection function, initial power spectrum, 
etc.. This might be a hard enterprise since it requires complex 
modellin g and exhaustive chec ks with numerical simulations as 
shown by lSanchez et al.l (l2008l) , all of them far beyond the scope 
of this paper. Nor do we intend to relate the measured features 
with scales that resemble particular events in the history of the 
universe, such as the sound h orizon at recombination as some 
authors have already done (e.g. lGuzik et al. 2007), since this can 
be misleading as pointed out by lSanchez et al.1 (l2008k . 

Our analysis consists of comparing the positions of features 
(i.e. peaks and valleys) in the linear regime of the quasar two- 
point correlation function, ^(r), at diff'erent redshifts. We are in- 
terested in a comparison among ^(r)s at diff'erent times over a 
wide range of scales, and not just finding an agreement on a par- 
ticular comoving scale as has been done in previous works (e.g. 
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Fig. 2. S/N analysis for the two redshift bin case. Upper panel: 
The percentage of ^(r) at large scales with S/N > 1, aver- 
aged over two redshift bins, i.e., the percentage of scales in 
^(r) we use to derive cosmological parameters. Lower panel: 
Averaged S /N map for the ^(r) in two redshift bins including 
only those scales where S/N > 1. Both panels corresponds to 
scales 100 h'^Mpc < r < 300 h'^Mpc and show higher values 
with lighter colours. 



R02). Regarding the number of redshift bins in which the quasar 
sample will be split, our approach is twofold: first, we divide the 
sample into two equal number redshift bins; second, we study 
the possible eff'ects of evolution in our results by dividing them 
into four redshift bins. The first choice aims to construct sam- 
ples of quasars as large as possible to allow for a fair detection 
of features in the ^(r) with signal-to-noise (S /N) ratios as high 
as possible for a wide range of comoving scales. Nevertheless, it 
should be noticed that the size of the redshift bins could allow for 



evolution of the quasar population within a single bin and thus 
bias the results. Therefore, we intend to understand possible evo- 
lution biases on the results that we found in the two redshift bin 
case by means of a four redshift bin analysis, even though this 
case will be aff'ected by larger errors bars in the ^(r). 

3.1. Computing the quasar correlation function 

In order to keep the noise as low as possible, we calculate the co- 
moving separations directly in three-dimensional, flat, comoving 
space, assuming a perfectly homogeneous metric. That is, we do 
not attempt the commonly used procedure that consists of deriv- 
ing the spatial two-point correlation function from the projected 
correlation function since that would require us to split pairs into 
a two-dimensional array, which would increase shot noise. 

The cosmology is implicitly involved in the spatial two- 
point correlation function through the computation of the co- 
moving distances. The comoving distance-redshift relation in a 
quintessence flat universe is given by: 



r(z) 



-f 

^0 Jo 



dz 



(1) 



where Dq = 1 - Dm and wq is constant. We explore diflTer- 
ent cosmologies, varying the parameters 0.1 < Qm ^ 0-55 and 
-1.1 < Wq < -0.1 in steps of ADm = 0.025 and Awq = 0.05. 
We are not interested in higher Qm values since there is com- 
pelling evide nce that we are living in a low m atter density uni- 
verse (e.g. iSanchez et al.ll2006l:[Spergel et al ]|1)07). 

For each pair of parameters, Qm and wq, we assign comov- 
ing positions of quasars by using their redshi fts in Eg. [D and 
compu te the correlation functions according to iLandy & Szalayl 



DD(r) - 2DR(r) RR(r) 
RR(r) ' 



(2) 



where DD(r), DR{r) and RR{r) are the normalised number of 
quasar-quasar, quasar-random and random-random pairs with 
comoving separations within r ± Ar/2. For the purposes of 
this work we use linear bins in comoving distances with Ar = 
5riMpc. 

The random sample was constructed to have the same angu- 
lar coverage of the survey and a smooth redshift distribution that 
matches that of the quasars. The angular mask was constructed 
using routines included in the HEALPi43 package to pixelize 
the sky coverage of the SDSS DR6. As regards the redshift dis- 
tribution of the random points, we have found that a remarkably 
good fit to the redshift distribution of quasars is the sum of five 
Gaussian functions as shown in the upper panel of Fig. [T] By 
fitting this function we have a smooth redshift distribution that 
closely follows that of the quasars. We then draw the redshifts 
of the individual random points from this distribution and obtain 
the distribution shown in the lower panel of Fig.[T] The random 
sample is 20 times denser than the sample of quasars in order to 
reduce shot noise. 

For the comoving scales that we are dealing with (r > 
100 /z"^Mpc), we compute ^(r) error bars using the bootstrap 
technique. This method becomes an unbiased estimate of the 
ensemble error when the number of pair coun ts is large com- 
pared to the numbers of objects in the sample (iMo et al 



^ H ierarchical Equal Area iso-Latitude Pixelization, iGorski et"aLl 
([2005h 
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Fig. 3. - wq) maps for the ''cross-correlation'' (CC^^^^) of the ^(r) in two redshift bins. The upper left panel shows CC^ 2] 
taking into account peaks and valleys in the ^(r), while the lower panels show the CC^^ only for peaks (left) and valleys (right). 
The cross-correlation in these panels were computed only for S/N > 1 , except for the upper right panel, where CC|^ is shown 
taking into account peaks and valleys without any S /N restriction. Higher values of the cross-correlation correspond to lighter 
colours. Labels in contour levels show the percentage of the map outside the contours. 



IShanks&Bovl3ll994 . Hence, for each correlation function, we 
also have 20 bootstrap estimates. 

In order to avoid spurious features, we smooth the correla- 
tion function and each bootstrap estimation using a Gaussian fil- 
ter with (Tsmooth = 5h~^Mpc. It is worth mentioning that smooth- 
ing the pair counts instead of the correlation function produces 
comparable results. We use these smoothed correlation functions 
to extract cosmological information and the smoothed bootstrap 
estimates to compute error bars. 

3.2. Estimating Dm and wq 

As i t has been stated in previous works (e.g. iRoukemal 
I2OOII R02), features at large comoving scales should remain un- 
changed over time for the proper cosmology. Our purpose is to 
extract as much information as possible from the features mea- 



sured in the ^(r)s at diff'erent redshifts. To achieve this, we do a 
full scale comparison among ^(r)s. This approach should intro- 
duce tighter constraints on the cosmological models. 

3.2.1. The statistical tool 

To quantify the goodness of each cosmological model pair 
(11m, >^q), we use a function that allows us to check for sim- 
ilarities among ^(r)s at diff'erent redshifts. This function is the 
cross correlation and it is an integral that expresses the amount 
of overlap between two functions. Formally, the cross correla- 
tion of two real functions / and g is defined as 



(f^8)(r)= I f(t)g(T^t)dt (3) 
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Fig. 4. (flu - wq) map shifted to lower values of wq. The map 
shows CC^^^^ of the ^(r) for two redshift bins. The minimum 
comoving distance of the map is rmin =110 h'^Mpc. Each model 
is computed only for S /N > 1. Labels in contour levels show the 
percentage of the map outside the contours. 

where r is a lag applied to g. Particularly, we are interested in 
the case of a lag r = 0. Hence, the discrete formula for the cross 
correlation in our study is 

CCiij] = (^i * ^j) = £ Ur) ^iir) (4) 

r-r,nin 

where ^ is the two-point correlation function and the subscripts 
/ and j denote particular redshift bins. It should be taken into 
account that we have restricted the sum over a finite range of 
comoving distances (r^/«, Tmax) at which we perform our analy- 
sis of the ^(r)s. Particularly, we choose as a maximum distance 
rmax = 300 h'^Mpc since larger scales have very low S /N, while 
rmin is used as a variable taking values > 100 h'^Mpc. 

Since we are interested in considering only the features that 
show a true signal, from now on (unless specified otherwise) we 
only use those points in ^(r) whose S /N ratio is greater than 
unity. Points that do not exceed this threshold are set to zero in 
the computation of the quantity CC. 

3.2.2. The two redshift bin case 

Our first approach is to split the quasar redshift distribution into 
two equal number subsamples. These subsamples are called the 
low redshift subsample, < z < 1.45, and the high redshift one, 
1.45 < z < 2.4 with median redshifts of 1.07 and 1.78 respec- 
tively. The purpose of this choice is to produce large samples of 
quasars that allow the measurement of many statistically signif- 
icant features. 

Which cosmologies are the most likely to produce signifi- 
cant features in ^(r) ? In order to answer this question we con- 
struct two significance maps which are shown in Fig. [21 The 
first map (upper panel) shows the percentage of scales in the 
^(r) that shows S /N > I (significant features) averaged over 
the two redshift bins, while the second one (lower panel) shows 
the averaged S /N for the ^(r) for two redshift bins, also using 



ing flat quintessential universes 5 

n,=0.175 w=-0.95 n,=0.250 w=-1.00 

M Q M Q 




100 200 100 200 

r (h^^ Mpc, comoving) 



Fig. 5. The two point correlation functions for the preferred (left 
panels) and standard (right panels) cosmological models. Upper 
and lower panels display the ^(r)s for the highest and lowest 
redshift bins respectively. Error bars were computed using the 
bootstrap re- sampling technique and are shown as a grey region. 
Up and down arrows point out features with S /N > I that remain 
fixed in comoving scales for both redshift bins. 



significant features only. The ^(r)s are restricted to the range 
100/z-^Mpc < r < 300/z-^Mpc. In both panels of Fig. [Slighter 
colours indicate higher values. From the first map we observe 
that our restriction on S /N implies that the useful range of scales 
on the ^(r), on average, runs from 23% to 33%. On the other 
hand, the second map shows that the averaged S /N lies in the 
range 1.39 to 1.51. From the contour levels and colours in this 
Fig. we can observe that there are many regions in both maps 
that are prone to a significant range of scales that are useful and 
a relatively high average S IN. Beyond that, an important issue 
that should be noticed is that all cosmological pairs on the map 
have useful scales in the ^(r) and therefore all models are plau- 
sible in showing coincident features at diff'erent times. 

To quantify the level of agreement among ^(r)s, we use the 
cross correlation denoted by CC^ where we have added an 
upper- script, (2z), denoting the current case. This statistic has 
the capability of enhancing features that occur at similar scales 
while erasing those that do not appear at the same distance. 
Moreover, since there is a sum over a wide range of comoving 
scales, several coincident signals for a particular cosmological 
pair produce higher values of CC implying a better match be- 
tween the ^(r)s at diff'erent redshifts. 

The CC^2] maps with Vmin = HO /z"^Mpc are shown in 
Fig. [3l All panels show contour levels around the CC^ 2] max- 
ima that indicate the percentage of the map outside their bound- 
aries. The upper left panel shows the CC^ 2] computed taking 
into account all features (peaks and valleys) in the ^(r)s. From 
this panel, we can see that the most likely cosmological models 
are in a small region in the lower left corner of the map. The 
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Fig. 6. S /N analysis for the four redshift bin case. Left panels: The percentage of the ^(r) at large scales with S /N > 1, averaged 
over the redshift bins 1-3 (upper panel) and 2-4 (lower panel). Right panels: Averaged S/N map for the f(r) for redshift bins 
1 and 3 (upper panel) and 2 and 4 (lower panel) including only those scales where S/N > 1. All panels correspond to scales 
100 h'^Mpc < r < 300 h'^Mpc and show higher values with lighter colours. 



black cross shows the maximum CC^^ value that corresponds 
to = 0.175 and wq = -0.95. 

Since this result is achieved by using all features in the ^(r)s, 
an interesting question is whether a particular feature (peaks or 
valleys) contributes to a particular region in the map. This issue 
can be answered by computing the cross correlation statistics 
using eq. (H but taking into account positive (peaks) or negative 
(valleys) significant (S /N > 1) values of the ^(r) only. The re- 
sulting maps can be seen on the lower left (peaks only) and lower 
right (valleys only) panels in Fig. [3] 

From the left panel, we can see that peaks restrict the area 
displayed by the contour levels compared with the ones obtained 
with all features. The contours are reduced, ruling out some cos- 
mological pairs with wq > -0.6 and are very slightly shifted 
to the right (higher values of Qm). Even so, the preferred cos- 
mological pair remains unchanged. On the other hand, valleys 
alone give lower flu values (~ 0.150) and produce an elonga- 
tion of contour levels towards higher values of wq. Nevertheless, 



their contribution to the combined feature analysis is not strong 
enough to prevail over the signal preferred by peaks. 

Finally, it is interesting to know what would be the eff'ect of 
relaxing the restriction imposed on the S /N of the features. In the 
upper right panel of Fig. [3] we construct the CC^ without 
imposing any value on the S/N of the ^(r)s. From the colour 
intensity of the map it is clear that the main eff'ect (compared 
to the upper left panel) is the widening of the area of preferred 
cosmologies. Even when this map enhances the importance of 
some cosmological models that are not statistically significant, 
the resulting preferred cosmology turns out to be very similar to 
the one obtained from significant features only. 

Since our result on wq is close to the limit of the wq range 
used so far, it would be interesting to test whether the exclusion 
of Wq < -1.1 is biasing our results. Therefore, we repeated the 
cross correlation analysis over a shifted range, -1.45 < wq < 
-0.50. As can be seen from Fig. [4] our results are stable, hence 
no bias arises from our original choice for the wq range. 
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To illustrate the ability of our method to compare ^(r)s at dif- 
ferent times, in Fig. [5] we show a comparison between the best 
cosmological pair obtained using the CC^^^j (^m = 0.175 
and wq = -0.95, left panels) and the standard model (AC DM, 
Dm = 0.250 and wq = -1.00, right panels). The grey region in 
each panel shows the Icr bootstrap errors bars. The best model 
is clearly better than the standard one given the number of coin- 
cident significant features at diff'erent scales. The arrows in this 
Fig. show features observed at the same comoving scale that sat- 
isfy S /N > 1 . While the preferred cosmological model obtained 
from our maps has four coincident features, we can only observe 
one matching feature for the standard cosmological model. 

We have also checked whether the resulting best model de- 
pends on the choice of rmin by computing CC^^ maps varying 
rmin from 120 to 190 h'^Mpc. In all cases we found that the pre- 
ferred pair of cosmological parameters is the same obtained for 
rmin = 110 h~^Mpc, indicating the stability of our results. 

3.2.3. The four redshift bin case 

Given the size of the redshift bins used in the previous case, 
the evolution of the quasar population within each bin might be 
biasing the results. In this subsection, we use four redshift bins 
to avoid too much evolution within a single bin. This analysis 
tests the stability of the results obtained with two redshift bins 
even when narrower redshift bins imply larger errors in the ^(r). 
We divide the sample into four equal number redshift bins: 

1. 0.00 <z< 1.07 

2. 1.07 <z< 1.45 

3. 1.45 <z< 1.78 

4. 1.78 <z< 2.40. 

With these subsamples we choose to calculate two estimates 
of the cosmological parameters using the cross correlations be- 
tween redshift bins 1-3 and 2-4. This choice will give us two es- 
timates statistically independent of each other, and, at the same 
time, by interleaving redshift bins, minimise evolutionary ef- 
fects. This approach will also allow us to estimate random errors 
in the resulting cosmological parameters by computing the stan- 
dard error in the mean between the two estimates. The resulting 
standard error can also be used as a fair estimate of random er- 
rors in the two redshift bin case. 

After computing the functions ^(r), we show in Fig. [6] the 
significance maps for the redshift bins 1-3 (upper panels) and 2- 
4 (lower panels). The percentage maps (left column) are shown 
for significant features averaged over two redshift bins. When 
combining redshift bins 1-3 (upper panel) we observe that 18%- 
26% is the range of scales for the ^(r)'s that are useful for our 
study while 19%-31% is the range of scales obtained for redshift 
bins 2-4 (lower panel). On the other hand, when analysing the 
averaged S /N maps (right column) for significant features only, 
the redshift bins 1-3 (upper panel) show a variation from 1.25 
to 1.59 while the average over redshifts 2-4 (lower panel) lies in 
the range 1.46 to 1.56. It is clear that all cosmological pairs have 
some portion of the ^(r) useful for our study. 

These S /N maps span similar value ranges to those corre- 
sponding to the two redshift bin case (see Fig. O. Even when 
computing the ^(r)s within narrower redshift bins produces 
larger error bars, on average, the S /N show no degradation com- 
pared to the previous case since some peaks and valleys in ^(r) 
are larger now. This could be a consequence of less evolution 
within a single redshift bin. Under this assumption, some fea- 
tures in the ^(r) in the two redshift bin case could have been 
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Fig. 7. Comparison of the ^(r)s between the two (blue and red 
curves) and four (grey areas) redshift bin cases for the preferred 
model obtained in the two redshift bin case. Error bars were 
computed using the bootstrap re- sampling technique. Each panel 
quotes the corresponding redshift ranges. 

eroded away by evolution. In Fig. [7] we observe a comparison 
between the ^(r)s obtained in the two (blue and red curves) and 
four (grey areas) redshift bins cases for the preferred cosmolog- 
ical model found in the previous subsection. Examples of our 
statement are the regions with 200 h'^Mpc < r < 250 h'^Mpc 
in the two upper panels, and the peak at ~ 130 h'^Mpc in the 
lower left panel of this Fig. 

Figure [8] shows the cross correlation maps corresponding to 
the calculation of equation |4] for redshift bins 1-3 (CC^^3p left 
panel) and 2-4 (CC|2 4], right panel). Each map is computed us- 
ing rmin =110 h'^Mpc. The preferred model in the CC[j^3^ map is 
Qm = 0.225 and wq = -0.90. On the other hand, the CC^^^l^ map 
better fits with Qm = 0.200 and wq = -0.95. We have checked 
that these results do not change when varying rmin from 120 to 
190 h'^Mpc. From these results the mean cosmological values 
and their corresponding standard errors are Qm = 0.21 ± 0.02 
and Wq = -0.93 ± 0.04. We observe that this result is consistent 
within a 1 (T interval with the one obtained for the two redshift 
bin case. 

In the light of our results, we consider as the most likely 
cosmological model to be the one obtained for the four redshift 
bin case since it is less aff'ected by evolution eff'ects while having 
similar S /N values to those obtained when using larger redshift 
bins. 

4. Summary and discussion 

In this work we have studied the evolution of the quasar two- 
point correlation function, ^(r), to constrain cosmological pa- 
rameters for flat quintessential universes, namely the matter adi- 
mensional density parameter Qm and the quintessence equation 
of state, Wq. For this purpose, we have drawn a large homo- 
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Fig. 8. (Qm - >^q) maps for the CC^"^^^ of the ^(r) for four redshift bins. The left panel show the CC^^^^^y while the right panels show 

CC|2 4]. The minimum comoving distance for each map is rmin =110 h'^Mpc. Each model is computed only for S /N > 1. Higher 
values of the cross-correlation correspond to lighter colours. Labels in contour level show the percentage of the map outside the 
contours. 



geneous quasar sample from the SDSS DR6, the largest quasar 
catalogue available at present. The size and homogeneity of the 
sample allow for highly reliable statisti cal studies of quasars . 
As in previo us studies of this aspect (e.g. Vallin otto et afll2007l: 
ISanchez etal. 2008: Croc ce & Scoccimarro 200 8|), we have fo- 
cused on analysing the position of features in ^(r), at large co- 
moving scales (r > 100/z"^Mpc) well into the linear regime of 
structure formation. The key point of the analysis lies in the as- 
sumption that, for the proper cosmological model, the comoving 
length scale of such features in ^(r) should remain fixed at dif- 
ferent stages in the evolution of the universe. In this work, we do 
not assume any particular scale at which the features (peaks and 
valleys) should be found, but we only require consistency among 
diff'erent redshifts ranges, imposing purely geometric constraints 
on the cosmological parameters. 

The study is focused on two diff'erent complementary ap- 
proaches: splitting the quasar redshift distribution into two and 
four redshifts bins. The first analysis is performed in order to 
have reliable estimations of ^(r) at large scales (where the signal 
is expected to be low) while the second one is performed aiming 
to rule out a possible bias in the previous results due to an evolu- 
tion of the quasar population within a particular redshift bin. The 
comparison among ^(r)s of diff'erent redshift bins is carried out 
by means of the cross-correlation function (CC). This particular 
statistic has the advantage of quantifying the amount of over- 
lap between two functions. In our case, the CCs are computed 
for modified ^(r)s, i.e, we have imposed a significance criterion 
which keeps the ^(r)s values when they satisfy the constraint 
\^(r)/cr(r)\ > 1, otherwise are set to zero. 

With regard to the two redshift bin case, the most likely 
model is Dm = 0.18 ± 0.02 and wq = -0.95 ± 0.04. The co- 
incidence among the results for the diff'erent reference scales 
(rmin) analysed is remarkable, implying that ^(r)s for both red- 
shift bins agree over a wide range of comoving scales, from 
~ 100 to ~ 300/z"^Mpc. From the four redshift bin case we ob- 
tain that the preferred cosmological model is Qm = 0.21 ± 0.02 
and Wq = -0.93 ± 0.04. When analysing the four redshift bin 
case, it becomes clear that the two redshift bin analysis is af- 



fected by evolution. The main observable eff'ect of evolution is 
that, on average, the amplitude of some features in the f (r) has 
been diminished, and as a consequence, the best model has been 
slightly shifted towards lower values of flu- 

The closest reference in the literature of an analysis similar to 
that performed in this work is the study of R02. Some of the sim- 
ilarities with this work are: we both only use a sample of quasars 
to constrain cosmological parameters, and neither of us search 
for a particular comoving scale. Beyond quasar sample sizes, 
the methodology diff'erences can be enumerated as follows: 

- We extract information from the position of peaks and val- 
leys in the ^(r)s while R02 used peak information only. 

- We restrict the ^(r)s with a significance criterion in order to 
avoid spurious detections while R02 used all the information 
without any restriction on S /N ratio. 

- The statistical tool proposed in this work intends to perform 
a full scale comparison among ^(r)s, i.e, to take into account 
more than one feature agreement among diff'erent redshift 
bins. On the other hand, R02 developed a method based on 
searching for the largest probability that the first local maxi- 
mum in ^(r) above a given reference scale occurs at the same 
comoving distance for consecutive redshift bins. 

- By construction, our method gives larger weights to stronger 
features while the probability method of R02 only considers 
the comoving length scale at which a peak is present. 

- We performed an analysis involving two and four redshift 
bins and found that evolution cannot be neglected. R02 per- 
formed their analysis using three redshift bins, extracting 
information from consecutive redshift bins. Their approach 
can still be aff'ected by evolution eff'ects. 

In their work, R02 found a matter density parameter in good 
agreement with our findings (see Table [T]). Nevertheless, they 
cannot impose strong constraints on the quintessence parameter, 
finding only an upper limit (wq < -0.50). Therefore, our results 
show that using the full large scale shape of the ^(r)s provides 
tighter constraints on cosmological parameters rather than bas- 
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Table 1. Constraints on cosmological parameters Qm and wq found in the literature under the assumption of perfectly flat universes. 
The last column quotes the data sets used in each work. 



Author 



Data 



lOK 2dFQS0^ 
WMAP2 + SDSS MGS^ 
WMAP + SDSS MGS + LRG^ 
RO45 WMAP36 + SNIa + LRG 
A06^ WMAP3 + SNIa + LRG 
WMAP + 2dFGRS^ 
WMAP3 + SNLS'^ 
DR6 SDSS QSO 



Roukema etal. (2002) 
Xesmark et al. (2004) 
^isenstein et al. (2005) 
Wang & Mukheriee (2006) 
Wang & M ukheriee (2006) 
Sanchez et a l. (2006) 
Spergel et al. (2007) 
This work 



0.25 + 0.10 
0.329 + 0.074 
0.326 + 0.037 



0.241 + 0.024 
0.238 + 0.015 
0.21 + 0.02 



-1.0,-0.5] 
-0.92 + 0.30 
-0.80 + 0.18 

-0 +0.206 
^•OO-^_0.iii -0.227 
_n QQQ+0-082 +0.159 
-rQfPSS -0.168 

-0.85!^]^ 
-0.967 + 0.073 
-0.93 + 0.04 



1: Quasars from the IQK release of the 2 degree Field Quasar Redshift Survev JCroom et 

2: First year Wilkinson Microwave Anisotropy Probe. ,.Spergel et al.. 1 20031) . 

3: Main Galaxy Sample, Strauss et al. 1 2002). 

4: Luminous Red Galaxies, Eisenstein et al. ( 2 0Qlh . 

5: lRiess et all (2004). 

6: Third year Wilkinson Microwave Anisotropy Probe. fspergel et all i2007l) . 
7: lAstieretaD [20061) . 

8: 2 degree Field Galaxy Redshift Survev. lColless et al.l 120011) . 
9: Supemovae Legacy Survey, lAstier et al.1 i2006l) . 



ing the analysis on the coincidence of a single particular scale 
only. 

From the comparison with results from a variety of data 
and techniques (see Table [TJ, we obs erve that our estimatio n of 
Qm is in good agre ement with R02, ISanchez et al.l (l2Q06h and 
ISpergel et al.l (l2007l) . In addition, our result for wq agrees with 
all sources quoted in Table [T] and has one of the smallest uncer- 
tainty intervals. 

What are the possible implications of our results in terms 
of structure formation compared to the initial conditions estab- 
lished by the ACDM concordance model? We observe a small 
discrepancy in the mean values obtained for the Qm and wq 
compared with a standard universe with cosmological constant 
(Qm = 0.25 and wq = -1.0). Our results are prone to prefer 
cosmologies with slightly lower values of Qm and higher values 
of Wq than those of the concordance model. Both tendencies ob- 
served in the parameters have similar consequences for structure 
formation. If all cosmological models were normalised to pre- 
dict the observed amplitude of density fluctuations at present, 
lower values of density woi ild predict hig her levels of density 
fluctuations in the past (e.g. ISpergeilll998l) . Under the same as- 
sumptions of density perturbation normalisation, earlier cloud 
formation and higher core densities are observed if the dynam- 
ics of the dark energy is enhanced, i.e. for m odels represented 
by an equation of state p arameter w n > - 1 (iBartelmann et al.l 
l2QQa iDolag et al.l [20041 : iMaio et al.l l2006h . Hence, at a given 
mass and redshift our results would predict higher abundance 
of more concentrated halos than in the model with the cosmo- 
logical constant. 

Our results rely on the analysis of quasar clustering alone, 
demonstrating that it is plausible to put strong constraints on 
cosmological parameters with only one suitable data set and a 
proper statistical treatment. Since our results are independent 
of the cosmic microwave background, type la supernovae and 
galaxy data, then large quasar samples in combination with all 
these data sets should imply a better determination of Qm and 

Wq. 



Acknowledgements 

We thank to the anonymous referee for important sugges- 
tions that have greatly improved the original manuscript. 
We also thank E.D. and Diego G. Lambas for carefully 
reading the manuscript and suggestions. HJM acknowledges 
the support of a Young Researchers' grant from Agenda 
Nacional de Promocion Cientifica y Tecnologica Argentina, 
PICT 2005/38087. This work has been partially supported by 
grants from Consejo de Investigaciones Cientificas y Tecnicas 
de la Republica Argentina (CONICET). Funding for the Sloan 
Digital Sky Survey (SDSS) has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the National 
Aeronautics and Space Administration, the National Science 
Foundation, the U.S. Department of Energy, the Japanese 
Monbukagakusho, and the Max Planck Society. The SDSS 
Web site is http://www.sdss.org/, The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Participating 
Institutions. The Participating Institutions are The University 
of Chicago, Fermilab, the Institute for Advanced Study, the 
Japan Participation Group, The Johns Hopkins University, the 
Korean Scientist Group, Los Alamos National Laboratory, 
the Max Planck Institut fiir Astronomic (MPIA), the Max 
Planck Institut fiir Astrophysik (MPA), New Mexico State 
University, University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observatory, and 
the University of Washington. 



References 

Adelman-McCarthy, J. K., et al. 2008, ApJS, 175, 297 
Astier, P., et al. 2006, A&A, 447, 31 

Bartelmann, M., Perrotta, R, & Baccigalupi, C. 2002, A&A, 396, 21 

Blanton, M. R., Lin, H., Lupton, R. H., Maley, F. M., Young, N., Zehavi, I., & 

Loveday, J. 2003, AJ, 125, 2276 
Caldwell, R. R., Dave, R., & Steinhardt, R J. 1998, Physical Review Letters, 80, 

1582 

Chen, G., & Ratra, B. 2004, ApJ, 612, LI 
CoUess, M., et al. 2001, MNRAS, 328, 1039 
Crocce, M., & Scoccimarro, R. 2008, Phys. Rev. D, 77, 023533 
Croom, S. M., Smith, R. J., Boyle, B. J., Shanks, T., Loaring, N. S., Miller, L., 
& Lewis, I. J. 2001, MNRAS, 322, L29 



10 



Zandivarez & Martmez: Constraining flat quintessential 



Daly, R. A., Mory, M. P., O'Dea, C. P., Kharb , P., B aum, S., Guerra, E. J., & 
Djorgovski, S. G. 2007, ArXiv e-prints, 710, arXiv:0 710.5112 

Dolag, K., Bartelmann, M., Perrotta, R, Baccigalupi, C, Moscardini, L., 
Meneghetti, M., & Tormen, G. 2004, A&A, 416, 853 

Dore, O., et al. 2007, ArXiv e-prints, 712, arXiv:0712. l399] 

Eisenstein, D. J., et al. 2001, AJ, 122, 2267 

Eisenstein, D. J., et al. 2005, ApJ, 633, 560 

Gong, Y., Wu, Q., & Wang, A. 2008, ApJ, 681, 27. 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, F. K., Reinecke, 

M., & Bartelmann, M. 2005, ApJ, 622, 759 
Gunn, J. E., et al. 1998, AJ, 116, 3040 

Guzik, J., Bernstein, G., & Smith, R. E. 2007, MNRAS, 375, 1329 
Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 

La Vacca, G., & Colombo, L. P. L. 2008, Journal of Cosmology and Astro- 
Particle Physics, 4, 7 

Lima, J. A. S., Jesus, J. R, & Cunha, J. V. 2007, ArXiv e-prints, 709, 
arXiv:0709.2195 

Maio, U., Dolag, K., Meneghetti, M., Moscardini, L., Yoshida, N., Baccigalupi, 

C, Bartelmann, M., & Perrotta, F. 2006, MNRAS, 373, 869 
Mo, H. J., Jing, Y. R, & Borner, G. 1992, ApJ, 392, 452. 
Peebles, R J. E. 1984, ApJ, 284, 439 
Perlmutter, S., et al. 1999, ApJ, 517, 565 
Ratra, B., & Peebles, R J. E. 1988, Phys. Rev. D, 37, 3406 
Richards, G. T., et al. 2002, AJ, 123, 2945 
Riess, A. G., et al. 1998, AJ, 116, 1009 
Riess, A. G., et al. 2004, ApJ, 607, 665 
Roukema, B. R 2001, A&A, 369, 729 

Roukema, B. R, Mamon, G. A., & Bajtlik, S. 2002, A&A, 382, 397 

Sanchez, A. G., Baugh, C. M., Percival, W. J., Peacock, J. A., Padilla, N. D., 

Cole, S., Frenk, C. S., & Norberg, R 2006, MNRAS, 366, 189 
Sanchez, A. G., Baugh, C. M., & Angulo, R. 2008, ArXiv e-prints, 804, 

arXiv:0804.0233 

Sapone, D., & Amendola, L. 2007, ArXiv e-prints, 709, arXiv:070 93792] 
Santos, R. C, & Lima, J. A. S. 2008, Phys. Rev. D, 77, 083505 
Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 
Schneider, D. R, et al. 2005, AJ, 130, 367 
Sen, A. A. 2008, Phys. Rev. D, 77, 043508 
Shanks, T., & Boyle, B. J. 1994, MNRAS, 271, 753. 
Spergel, D. N. 1998, Classical and Quantum Gravity, 15, 2589 
Spergel, D. N., et al. 2003, ApJS, 148, 175 
Spergel, D. N., et al. 2007, ApJS, 170, 377 
Strauss, M. A., et al. 2002, AJ, 124, 1810 

Szydlowski, M., Hrycyna, O., & Kurek, A. 2008, Phys. Rev. D, 77, 027302 
Tegmark, M., et al. 2004, Phys. Rev. D, 69, 103501 

Vallinotto, A., Dodelson, S., Schimd, C, & Uzan, J.-R 2007, Phys. Rev. D, 75, 
103509 

Wang, Y, & Mukherjee, R 2006, ApJ, 650, 1 

Wang, Y, & Mukherjee, R 2007, Phys. Rev. D, 76, 103533 

Wei, H., & Zhang, S. N. 2008, ArXiv e-prints, 803 JarXiv:0803.3292l 
Yahata, K., et al. 2005, PASJ, 57, 529 
York, D. G., et al. 2000, AJ, 120, 1579 

Zhang, Q.-J., Cheng, L.-M., & Wu, Y.-L. 2007, ArXiv e-prints, 708, 

arXiv:0708.2164 

Zhu, Z.-H., & Sereno, M. 2008, ArXiv e-prints, 804 JarXiv:0804.2917l 



